function y = Cn_deltarr(x)
 
 if x >= 0
     c = [0 -6.7535 6.8214 -2.2465 0.2953 -0.0418 1E-12];
     p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
 else
     c = [0 -6.7535 -6.8214 -2.2465 -0.2953 -0.0418 -1E-13];
     p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
 end
 
  y = c * p;
  
end
